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This paper presents the nonlinear dynamics of laser cooled 
and trapped cesium atoms placed inside an optical cavity and 
interacting with a probe light beam slightly detuned from the 
6S1/2F = 4 to 6P3/2F = 5 transition. The system exhibits 
very strong bistability and instabilities. The origin of the 
latter is found to be a competition between optical pumping 
and non-linearities due to saturation of the optical transition. 

PACS: 32.80.Bx, 32.80.Pj, 42.65.Pc 



I. INTRODUCTION 

In an atomic vapor, the interaction of the atoms with 
one or several near resonant electromagnetic fields is com- 
plicated by the fact that the various velocity classes have 
different detunings from the fields, except when this de- 
tuning is very large. Laser cooled atoms in a magneto- 
optic trap jjj can have velocities as low as a centime- 
ter per second. This means that their Doppler width is 
smaller than the natural linewidth and that a laser field 
can be set close to resonance with an atomic transition, 
all the atoms having the same detuning from the field. 
In such conditions, the interaction between the atoms 
and the field is well characterized and one can take ad- 
vantage of the strong non-linearities of atomic systems 
while keeping the absorption rather small. In particular, 
it was shown that a probe beam going through a cloud 
of cold atoms could experience a strong gain due to Ra- 
man transitions involving the trapping beams Q . When 
the atoms are placed in a resonant optical cavity without 
a probe beam, laser action corresponding to that gain 
feature was demonstrated M. 

When cold atoms interact with a probe laser beam 
inside an optical cavity, bistability is easily observed at 
very low input powers (as low as 5/iW) B. This bista- 
bility effect was observed in the presence of the cooling 
beams. However, to investigate the nonlinear dynamics 
of a collection of cold atoms in a cavity in more detail, 
we needed better controled conditions, and we studied 
the behaviour of the system in the absence of cooling 
beams, right after the trap is turned off. In that case, in 
addition to bistability, new features were found, such as 
very pronounced self pulsing oscillations in a wide range 
of experimental parameters. 



While instabilities have been observed in similar condi- 
tions in atomic vapors the present situation is much 
easier to analyse, and we show hereafter that we have 
been able to find a rather simple model for these instabili- 
ties, that invokes a competition between optical pumping 
and saturation. 

The exprimental procedure is described in section 2. 
In section 3 we give a model explaining the dynamical 
behaviour of the system and compare its results with the 
measurements. 



II. BISTABILITY AND INSTABILITIES WITH 
COLD ATOMS 

In the experiments described in the following, we pre- 
pare a cloud of cold cesium atoms in a cell using the 
background pressure to fill the trap. The trap operates in 
the standard way Q , with three orthogonal circularly po- 
larized trapping beams generated by a TkSapphire laser 
and an inhomogeneous magnetic field. The Ti:Sapphirc 
laser is detuned by 2.5 T (T being the linewidth of the 
upper state) on the low frequency side of the 6Si/ 2 F = 4 
to 6P3/2-F 1 = 5 transition. We obtain a cloud of cesium 
atoms the typical temperature of which is of the order of 
1 mK, which gives a Doppler width much smaller than 
the natural width. The diameter (2.5 cm) and power 
(20 mW/cm 2 ) of our trapping beams allow us to obtain 
large clouds (about 5 mm in diameter) with densities of 
the order of 10 10 atoms/cm 3 . The relevant parameter in 
the experiment is actually the number of trapped atoms 
in the probe beam, which is measured from the change in 
the intensity of the probe beam transmitted through the 
cavity with and without trapped atoms. This number 
is found to be ranging between 10 7 and 10 8 depending 
on the pressure of the background gas. As usual, the 
atoms non resonantly excited into the 6P3/2-F = 4 state 
and falling back into the F = 3 ground state are re- 
pumped into the cooling cycle by a laser diode tuned to 



the 6S 1/2 F 



3 to 6P 3/2 F 



4 transition. 
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The cavity is a 25 cm long linear asymetrical cavity, 
close to half-confocal, with a waist of 260/im. Because 
the cell has optical quality antireflecting windows, we 
can build a good finesse optical cavity around the atomic 
cloud (Fig.l). Losses due to the two windows are of the 
order of 1%. The input mirror has a transmission coeffi- 
cient of 10%, the end mirror is a highly reflecting mirror. 
The cavity is in the symmetry plane of the trap, making 
a 45° angle with the two trapping beams that propagate 
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in this plane. 
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FIG. 1. Experimental set-up showing the cell containing 
the cold atom cloud in an optical cavity; PBS: polarizing 
beamsplitter, QWP: quarter wave plate, PD1, PD2: photo- 
diodes; PD1 and PD2 measure the powers respectively trans- 
mitted and reflected by the cavity. 

To look for bistable behaviour of the optical cavity con- 
taining cold atoms, we send a circularly polarized probe 
beam into the cavity. It can be detuned by to 130 MHz 
on either side of the 6Si/ 2 F = 4 to &P^/ 2 F — 5 atomic 
transition. 

We measure the power of the beam transmitted 
through the cavity while scanning the cavity length for 
a fixed value of the input intensity, as shown in Fig. 2, or 
the input intensity for a fixed value of the detuning. The 
recording shows the characteristic hysteresis cycle due to 
bistability, where the output power switches abruptly be- 
tween low and high values when the length of the cavity 
is scanned. Switching and hysteresis were also observed 
when the input power is scanned at fixed cavity length. 
In some cases overshoot and oscillations in the output 
power were recorded B. 
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FIG. 2. Recording showing the bistable switching from 
low to high transmission and back when the cavity length 
is scanned across the cavity resonance. The trapping laser 
beams are on. The laser is detuned by 22r on the high 
frequency side of the atomic transition, the input power is 
100 fiW. 

The shapes of the curves in such a case can be grossly 
interpreted from the theory of bistability with two-level 
atoms. However, the fit is approximate and the the- 
ory, even including single mode instabilities , cannot 
explain the overshoot and oscillations ||. Actually the 
transition under investigation in cesium is far from being 



a two-level one in the presence of the trap. The fact that 
standard bistability is observed in most cases to be in fair 
agreement with two-level atom theory can be explained 
by the fact that the trapping beams randomize the the 
ground state population among the various Zeeman sub- 
levels. 

To thoroughly investigate the phenomenon, it was thus 
indispensable to study the system without the trapping 
beams. After the trap is loaded, we turn off the trapping 
laser beams in order to get unperturbed atoms. We have 
about 20 ms to perform the measurements before most of 
the atoms have escaped out of the interaction region due 
to free fall and expansion of the cloud. In the experiment, 
the bistability parameter C (see definition in section 3) 
can be as high as 400 just after the atoms have been 
released. 

In a broad range of experimental parameters, we ob- 
served bistability and instabilities that exhibit unusal fea- 
tures. Instabilities are present within the whole range of 
accessible detunings (from — 25r to 25r, where T is the 
linewidth of the excited state, r/27r = 5.2 MHz), and for 
input powers ranging from 50 to 300//W. These oscilla- 
tions are somewhat similar to the ones observed in the 
presence of trapping beams but those were obtained only 
at high intensity or small atomic detuning. Here, on the 
contrary, they are observed very easily. 

Fig. 3 shows a set of recordings of the output power of 
the cavity when the length is scanned, for different values 
of the input power. At very low power (not shown), the 
curves exhibit neither bistability nor instabilities. Start- 
ing for an input power of the order of 30/xW, oscillations 
appear on the left hand side of the cavity resonance curve, 
within some range of cavity detuning. This side is the one 
on which a bistable switching would occur in a saturated 
two-level atomic system. At higher input powers, the os- 
cillations disappear and only bistability persists. Let us 
mention that when the power of the probe beam is too 
high, the atoms are expelled very rapidly from the beam 
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FIG. 3. Recording of the output power P ou tof the cavity 
containing cold atoms when the cavity length is scanned, for 
four different values of the input power Pi n . The trapping 
beams are off. The atomic detunig is the same as in Fig. 2. 
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and only the very early part of the signal can be con- 
sidered as significant of the nonlinear behaviour of the 
atoms. 

Once enlarged these oscillations look clearly like self- 
pulsing. Their frequency is comprised between 100 KHz 
and a few MHz, and they are not due to the scan of the 
cavity length. To investigate them in more detail, one 
would want to record them at a fixed cavity length. How- 
ever, it is not possible to keep the optical cavity length 
perfectly constant in time because the atoms escape from 
the original cloud, thus changing the index of refraction. 
But in such conditions, the effective cavity scan is slow 
and the oscillations are observed on longer time dura- 
tions. Such a recording is shown in Fig. 4. 
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FIG. 4. Recording of the instabilities for a fixed geometri- 
cal cavity length. The optical length is slowly scanned by the 
decrease of the number of atoms. The detuning is the same 
as in Fig. 2, the input power is 80 /iW, the trapping beams 
are switched off at t = 0. 



III. MODEL FOR INSTABILITIES 

To understand these oscillations, one has to take into 
account the hyperfine and Zeeman structure of the con- 
sidered states. Various optical pumping effects can oc- 
cur and phenomena linked to it like bistability, multista- 
bility ggmi and instabilities §,§|||,[|^|^,|^,|f| 
have been predicted and observed in alkali vapors. Insta- 
bilities necessitate a strong coupling between the atoms 
and the field, that is a small detuning. Considering that 
the Doppler width is of the same order as the hyperfine 
structure in the ground and excited state of cesium, the 
system is rather intricate to describe, due to the compet- 
itive action of the various velocity classes and hyperfine 
sublevels. 

As already pointed out, the situation is much simpler 
with cold atoms, where one can consider that the field in- 
teracts with one hyperfine transition only. Owing to this, 
we have been able to develop a simple model to under- 
stand the origin of the observed oscillations. Roughly, as 
shown below, they result from the competition between a 



fast nonlinear process, the saturation of the optical tran- 
sition and a much slower one, the optical pumping. 

The saturation of the optical transition causes a de- 
crease of the linear index of refraction of the atomic 
medium when the intensity increases. On the contrary, 
optical pumping by circularly polarized light is at the 
origin of a non-linearity that increases the index of re- 
fraction with the light intensity: when the atoms are 
submitted to circularly polarized light, they tend to accu- 
mulate in the magnetic sublevels with high mp number, 
which have the largest coupling coefficient with the elec- 
tromagnetic field. Therefore, the two nonlinear processes 
have opposite effects and compete in our system. The re- 
laxation oscillations are a consequence of the significant 
difference in the characteristic times of theses processes. 

Optical pumping tends to empty out the magnetic 
sublevels of low magnetic number to accumulate all the 
atoms in the sublevels with highest magnetic number 
(to p = A to mp = 5 transition) . Due to the high number 
of magnetic substates, it takes a rather long time, start- 
ing from an equally distributed population, to complete 
the optical pumping to the highest mp sublevels of the 
ground and excited states. 

The time evolution of the sum of the populations of 
these two sublevels for several values of the Rabi fre- 
quency of the probe beam (calculated as an average over 
the hyperfine transitions) is shown in Fig. 5. One can see 
that even for values close to or larger than the saturation 
value, the optical pumping rate is much smaller than the 
natural linewidth. Thus the response time of the nonlin- 
ear susceptibility due to optical pumping is much larger 
than the one due to saturation of the optical transition, 
which, for a field detuned from resonance, can be consid- 
ered as being of the order of the detuning. 
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FIG. 5. Optical pumping as a function of time : the curves 
show the calculated variation of the sum N of the populations 
of the rriF = 4 sublevel of the ground state and of the mF~5 
sublevel of the excited state as a function of time, starting 
from a population equally distributed among the ground state 
sublevels, when the intensity I of the pumping light (defined 
by Eq. (1)) is equal to (from bottom to top) : 1, 5, 10, 12, 
15, 20, 30, 40, 60. Time is in units of F. The atomic detuning 
is 20F. 
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Although the behaviour of the system, involving many 
variables (all the hyperfine Zeeman populations and co- 
herences and the field), is quite complex, the underlying 
mechanism can be explained with a simple model, in- 
volving basically two differential equations that give the 
evolution of the intracavity field and that of the atomic 
orientation in the ground state. To write these equations, 
we first introduce our basic notations. 

The intracavity intensity I of the laser is normalized 
to the saturation intensity : 



g 2 H 2 
r 2 /4 



(i) 



where g is the coupling constant of the atoms with the 
field, 



d 2 ui L 
2e a frSc 



(2) 



d is the atomic dipole, lul the frequency of the probe 
laser and S the cross section area of the beam. \a\ 2 is 
the electric field squared, expressed in units of number of 
photons per second. 

The round trip phase shift <£> C av of the field in the cavity 
(assumed to be a ring cavity) is the sum of four contri- 
butions. First, the phase shift $o proportional to the 
geometrical length of the cavity. Second, we have two 
contributions due to the presence of atoms in the cavity, 
a linear phase shift, 



$ L = 2Ng 2 /T5 
and a nonlinear Kerr-like phase shift 



NL 



-KI. 



The nonlinear coefficient K is given by 
K = 4A^ 2 /r<5 3 , 



(3) 



(4) 



(5) 



where N is the number of atoms and S the detuning of the 
atomic transition frequency u>o from the field frequency 
wl, normalized to the atomic transition linewidth T/2 



5 = 2(c^o - WL)/r. 



(6) 



<E>l and $nl are the phase shifts corresponding to the 
presence of two level atoms in the cavity. If we now con- 
sider that the ground and excited states have several Zee- 
man sublevels, the main additional contribution when the 
atoms interact with circularly polarized light is a term <f> p 
coming from the change in the populations of the ground 
state sublevels and proportional to the ground state ori- 
entation p (normalized to 1). Since the square of the 
Clebsch-Gordan coefficient of the mp—^ to rap = 5 tran- 
sition is about twice the mean square Clebsch-Gordan 
coefficient for the F=4 to F=5 transition, $ p is equal to 
$l for p = 1, that is when the ground state is completely 
pumped. Thus, we can write 



The total phase shift in the cavity is then: 



(7) 



(8) 



The ground state orientation p increases with the in- 
tracavity intensity I at rate (51 and decays at a rate 7 P 
due to magnetic precession in transverse fields and to 
transitions to other hyperfine sublevels (via non resonant 
transitions): 



dp /At = -j p p + /3I(1 -p). 



(9) 



The pumping rate coefficient (3 is computed from the 
calculation presented in Fig. 5 and the relaxation rate 
7 P is evaluated from the experimental parameters. 

The change of the intracavity field a on a round trip 
of time duration t is due to the driving field a; n enter- 
ing through the coupling mirror of transmission t, to the 
mirror decay coefficient 7 cav (with 7 cav = t 2 /2) and to 
the round trip phase shift $ cav : 



rda/dt = to in - (7cav - «$cav)a. 



(10) 



The Kerr-like non-linearity has been assumed to have 
an instantaneous response. In the absence of optical 
pumping this system is well known to become bistable 
when the intensity is larger than a threshold intensity 
Jbist given by 



^bist - 8 7c 2 av /3v / 3K 



(11) 



In the presence of optical pumping eqs. (8) and (9) 
have to be solved numerically They involve two different 
rates: the optical pumping rate f3I and the field evolution 
rate in the cavity Jchv/t (7 cav /27rr w 5 MHz) which is 
usually larger than the optical pumping rate, except at 
very high pump powers. 

These equations have been used to calculate the mo- 
tion of the system. In some range of initial conditions, 
oscillations and limit cycles are found in the intracavity 
intensity as well as in the output intensity. To compare 
these results with the experimental data, we have cal- 
culated the output intensity when $o, i.e. the cavity 
length, is slowly scanned for several values of the input 
intensity. The result is shown in Fig. 6. It can be seen 
that the curves reproduce the experimental recordings in 
a satisfactory way. 

In particular, they show oscillations setting in for inter- 
mediate powers. When the input power is strong enough, 
the oscillations disappear, due to the fact that the opti- 
cal pumping is fast enough to bring all the atoms in the 
sublevels of high magnetic number before the oscillations 
can start. The calculations show that in the unstable 
region, the ground state orientation may vary by quanti- 
ties as small as 1%. As a consequence, the period of the 
instabilities can be much smaller than the typical optical 
pumping time shown in Fig. 5. 
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FIG. 6. Calculated output power of the cavity containing 
cold atoms when the cavity length is scanned for four values 
of the input power chosen below and above the bistability 
threshold Pbist (-Pbist is in the absence of optical pumping, cf 
eq.(10)). 

This oscillatory behaviour can be understood by con- 
sidering that the optical pumping changes the optical 
length of the cavity. It can then scan the cavity length 
back and forth in the vicinity of the bistable regime, caus- 
ing periodic abrupt changes in the cavity transmission as 
shown in Fig. 7. A similar behaviour was observed in a 
different system with thermal effects Q| . 

Instabilities can also be found in the absence of re- 
laxation for the orientation when the cavity length is 
scanned. This occurs for high intensities and large num- 
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FIG. 7. Diagramme showing the self-pulsing of the bistable 
system due to optical pumping (for the case of $l an d 
$ p < 0); (a) the cavity is in the high transmission state and 
optical puming increases the orientation and decreases $l; 
(b) abrupt switch towards the low transmission state; (c) the 
cavity is in the low transmission state, depumping is dom- 
inant and the orientation decreases; (d) switch towards the 
high transmission state. 



bers of atoms. In such a case, in the vicinity of a reso- 
nance, the optical pumping takes over, brings the cavity 
first into resonance and then beyond the resonance. After 
this overshoot, the cavity length scan brings the cavity 
back into resonance, the optical pumping starts again 
and so on. Thus the orientation increases by steps, each 
time the light enters the cavity. 

The validity of the model was checked by further ex- 
periments. When the atoms are released from the trap, 
it is possible to optically pump them into the ttif = 4 
magnetic sublevel of the ground state with an additional 
circularly polarized beam parallel to the probe beam, but 
closer to the atomic resonance. This prepumping is done 
in the presence of a magnetic field directed along the 
cavity. In such conditions the instabilities disappear. If 
the magnetic field is absent, the orientation created by 
the pump field is destroyed by the Larmor precession in 
tranverse magnetic fields and the instabilities persist. 

A more complete treatment of the atomic non-linearity 
and of the optical pumping was also performed, where ab- 
sorption and saturation of the optical non-linearity were 
taken into account. In this case, the separation of linear 
and nonlinear phase shifts is no longer possible. Instead, 
one defines a total phase shift due to two-level atoms 



*i = 



2Ng 2 



5 + i 



(12) 



r i + s 2 + 2i 

or with the help of the bistability parameter C given by 



C = g 2 N/ lc&v T, 



2Cj cllv (S + i) 
l + 6 2 + 2I ' 



(13) 
(14) 



$i now includes a contribution due to the linear ab- 
sortion and dispersion of the atoms. The contributions 
$l and 'J'nl introduced below are simply the first two 
real terms of the expansion of formula (11) in powers of 
I. The total phase shift, including optical pumping now 
writes 



$cav = $o + *i(l+p)- 



(15) 



For consistency, the saturation of the optical pumping 
was also taken into account: 



dp/dt 



- lpP + P l + S 2 + 2I {1 - p) - 



(16) 



The evolution of the system was calculated again us- 
ing this more elaborate model. This yielded only minor 
changes in the results, the general behaviour predicted 
by the simple model for the instabilities remaining the 
same. 

Let us note that for high powers, the probe laser tends 
to push the atoms out of the beam, especially when its 
frequency is close to resonance. This phenomenon could 
give rise to bistability due to the change of the effective 
linear index of refraction of the medium with the number 
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of atoms in the interaction zone. A careful study of the [13] G. Giusfredi, S. Cecchi, E. Pietrella, P. Salieri and F.T. 
behaviour of the cold atoms in the probe beam has shown Arecchi, in "Instabilities and Chaos in Quantum Optics" 

that such mechanical effects were negligible under our eds: F.T. Arecchi and R.G. Harrison, Springer (1987). 

experimental conditions. [14] S. L. McCall, Appl. Phys. Lett. 32 (1978) 284. 



IV. CONCLUSION 



The recent development of the magneto-optic trap, 
which enables to get clouds of motionless atoms with a 
density comparable to that of the atomic beams is of 
great interest for nonlinear optics. In such traps, one can 
get large nondinearities by setting the driving field close 
to resonance without having much absorption. The in- 
stabilities described and interpreted in this paper show 
that it is possible to get new nonlinear phenomena in well 
characterized conditions. These experiments open the 
way to nonlinear and quantum optics using cold atoms 
in resonant cavities. 
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